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Abstract 



We study a relation between two integrability conditions, namely the Yang-Baxter and 
the pair propagation equations, in 2D lattice models. While the two are equivalent in the 
8- vertex models, discrepancies appear in the 16- vertex models. As explicit examples, we find 
the exactly solvable 16-vertex models which do not satisfy the Yang-Baxter equations. 
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1. Introduction 

In last two decades, much progress has been made in 2D integrable systems both in 
lattice statistical models and in continuum field theories. Recently, this progress has been 
associated with beautiful mathematical structures such as universal Grassmann manifold [T| , 
Kac- Moody algebra p| and quantum group [||]. 

In 2D lattice models, there is one approach, which is based on transfer matrices (TMs) 
and it has been proved most successful. As Baxter showed, one can construct infinite num- 
ber of commuting conserved quantities through these TMs Q. A sufficient condition for the 
commutativity is that the Boltzmann weights of the 2D lattice models satisfy the famous 
Yang-Baxter equations (YBEs). There can obviously exist many exactly solvable models 
which do not satisfy YBEs. Since these are exactly solvable, one needs another scheme to 
solve these models if they exist. 

There is another approach, which is based on the so-called pair propagation equations 
(PPEs) appearing in the analysis of the algebraic Bethe ansatz. According to this method, 
the Boltzmann weights satisfy non-linear coupled equations. These equations become man- 
ageable if the Boltzmann weights are defined on some algebraic curves. 

In this paper, we want to study some 2D lattice models which can be exactly solvable 
while they do not satisfy YBEs. We are looking for our candidates from the 16-vertex 
models f||J. What we are going to show first is a relationship between YBEs and PPEs. 
Though YBEs and PPEs are equivalent in the 8-vertex model, discrepancies appear in 16- 
vertex models. Since YBEs restrict possible candidates so strongly, PPEs can cover more 
exactly solvable models which do not satisfy YBEs. We give explicit examples for which we 
compute exact eigenvalues of transfer matrices. 

2. The Pair Propagation and Conjugate Pair Propagation Equations 

We follow the notation of Baxter ^J]. The Boltzmann weights of the symmetric 16-vertex 
models are given by 

R(±, ±; ±, ±) = a, R(±, T; ±, T) = &, R(±, =f; t, ±) = c, R(±, ±; =F, t) = d, 

R{±, t; T, T) = e, R(±, ±; T, ±) = k, R(±, ±; ±, =f) = h, R(t, ±; T, T) = l (1) 

The Yang-Baxter equations are given in the forms; 

X R(p, C; v, P)R'(p, a; 0, R "(v, 0; v> <*) 
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= R"(»,p;v,<i>)R'((i>,C,v,l3)R(v,oc;v,0. (2) 

According to the Bethe ansatz, eigenf unctions y(Pi, p 2 , ■ ■-,Pn) of transfer matrices T(v) 
for N horizontal sites become in the forms of the direct products of each variables such as 
y(/3i, /3 2 , Pn) — 9i(Pi) ^92^2) <£> ••• ® 9n(Pn)- These are eigenf unctions of transfer matrix 
on the upper layer. We multiply these eigenf unctions to transfer matrices, and obtain 

(T(v)y) a = Tr (G 1 (a 1 )...G N (a N )) , with (G^a))^ = R (V, a] u, P)gi(P). (3) 



Explicit forms of Gj(±) are 
Gi{+) = 
Gi(-) = 



a-9i{+) + hgi(-) %(+) + %(-) 

egi(+) + cgi(-) bgi(+) + Zft(-) 

IgM + bgi(-) cg%(+) + egi(-) 

dgi{+) + kgi(-) %(+) + osse- 



in order to be solved exactly, it is necessary that there exist the a-independent pairs of 
matrices Pj,Pj+i, which transform Gi(a) into upper triangle forms; 



For simplicity, we choose det Pj = 1 for % = 1, 2, N and we parametrize them in the forms; 
Pj = ( , \ M \ V Then Eq. (HI) is written in the forms ; 

v **(-) / 

G i (a)P i+1 = P i H i (a), or Pr^a) = H^Pr^. (5) 
As Hi{a) are in upper triangle forms, we obtain; 

G '"»te-) h<°> ("<-!)• < 6) 

( -ft(-), Pi(+) ) Giipt) = g"{a) ( -p i+ i(-), p i+ i(+) ) . (7) 

We call Eq. as the pair propagation equations (I) and Eq. (|7|) as the conjugate pair 
propagation equations (I). By using P(//, a; u, P), the pair propagation equations (I) are 
given by 

^ l R(fjL,oi;v,P)g i (fi)pi + i(u) =^(a)pi(^). (8) 

0,v 
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Explicit forms of these equations are 

' (a</<(+) + hgi(-))p i+1 (+) + (k gi (+) + d gi (-))p i+l (-) = g' i (+)p i (+), 

(l 9i (+) + 6ft(-))p i+ i(+) + (cfc(+) + e ft (-))p m (-) = ft(-)p<(+), fg , 

(eft(+) + cgi(-))p i+1 (+) + (b 9i (+) + Z ft (-))p i+1 (-) = 1 1 
(%(+) + fc</i(-))p*+i(+) + (%(+) + o^(-))pi+i(-) = ^-(-)pi(-). 

While the conjugate pair propagation equations (I) are given by 

^2R(fi,a;u,/3)gi(/3)qi(fj.) = g- (a)q i+1 (v), (10) 

where we use notation q%(+) = —Pi(—), Qi{—) = Pi(+)- Explicit forms of these equations are 
obtained from Eq. © by replacing #-(±) -> g-(±), Pi+i(+) -> -Pi(-), Pi+i(~) Pi(+), 
Pi(+) Pi(-) Pi+i(+)> d, k. 

Next we consider the second type of the pair propagation equations, which we call the 
pair propagation equations (II). If the models are exactly solvable by using eigenfunctions 
y(Pi, @2, •••) Pn) — 5 , i(A)®5 , 2(/32)®---®5 , v(/9v) acting on upper layer of the transfer matrices, 
it is exacly solvable by using eigenfunction y(ax,ac2, ...,«») = <?i(«i) ® £b(«2) <8> ••• ®^v(«v) 
acting on lower layer. Similar equations corresponding to Eqs. (0),(||) are given by 

(yT{v)) p = Tr(G 1 {/3 1 )...G N {/3 N )) , with (g<(/3)) = «5 f, (11) 

and 

Pr^m+i = = ( *W ) • (12) 

The pair propagation equations (II) are given by 

^R(n,a]v,(3) gi (a)p i+l (v) = g'i(/3)pi(fi) . (13) 

By the symmetry of the Boltzmann weights, explicit forms of these pair propagation equa- 
tions (II) are obtained from the pair propagation equations (I) by replacing untilde variables 
into tilde variables and c d,h <-> I. Similarly we obtain the conjugate pair propagation 
equations (II) from Eq. ([!(]) by the same replacement. 

Our strategy to solve the pair propagation equations is the following. These pair propa- 
gation and conjugate pair propagation equations are special bilinear equations of four vari- 
ables such as <7i(±), (7j(±),pj(±),pj +1 (±), and it is difficult to solve these equations directly. 
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Then we first derive non-linear equations, where only two ratio of the variables such as 
r i — Pi{~ )/Pi(+)j r i+i — Pi+i(~ appear. Instead of solving four coupled equations 
in the pair propagation equations, we first solve these four non-linear equations. For each 
solution of these non-linear but two variable equations, only one of four equations of the pair 
propagation equations is independent, and from that we can obtain the eigenvalues of the 
transfer matrices. 

From condition that there exists non-trivial solutions for g>j(±), <?j(±) in Eq. we obtain 

rj + rj +1 - T^rl, + 1) - T 2 r t r l+1 + T 3 r t (l + rj +1 ) + r 4 r m (l + r 2 ) = (14) 
where 

T 1 = (cd - ek)/{ab - hi), T 2 = (a 2 + b 2 + e 2 + k 2 - c 2 - d 2 - h 2 - l 2 )/(ab - hi), 
T 3 = (cl + dh — ak — be) / (ab — hi) , Y± = (ae + bk — ch — dl) / (ab — hi) 

and 

n = Pi(-)/Pi(+), r i+1 = p i+ i(-)/p l+ i(+). 
From the condition that there exists non-trivial solutions for pi(±) , pi + \(±) , we obtain 

s 2 + s- - r 5 (s 2 s[ 2 + 1) - r 6SiS ; + r 7S ;(i + s 2 ) + r gSl (i + s f) = o (15) 

where 

T 5 = (cd - hl)/(ab - ek), T 6 = (a 2 + b 2 + h 2 + I 2 - c 2 - d 2 - e 2 - k 2 )/(ab - ek), 
Tj = (ce + dk — ah — bl) / (ab — ek) , Tg = (al + bh — ck — de) / (ab — ek) , 

and 

Si = gi(-)/gi(+), *1 = &(-)/&■(+)■ 

Simlarly, from the condition that there exists non-trivial solutions for ^(±),p i+1 (±),we 
obtain 

r 2 + s 2 - T 9 (r 2 s 2 + 1) - T 10 r iSl + r n n(l + s 2 ) + T 12Sl (l + rf) = (16) 
where 

T 9 = (bd - eh)/(ac - kl), T w = (a 2 + c 2 + e 2 + h 2 - b 2 - d 2 - k 2 - I 2 ) /(ac - kl), 
r n = (bl + dk — ah — ce) / (ac — kl), T 12 = (ae + ch — bk — dl) / (ac — kl). 

Finally, from the condition that there exists non-trivial solutions for <fc(±),pi(±), we 
obtain 

s? + r 2 i+1 - r 13 (sfr 2 + 1) - r 14 s-r m + r 15 s-(l + r 2 +1 ) + r 16 r m (l + sf ) = (17) 
where 

T 13 = (bd - kl)/(ac - eh), T 14 = (a 2 + c 2 + k 2 + I 2 - b 2 - d 2 - e 2 - h 2 )/(ac - eh), 
Ti5 = (be + dh — ak — cl) / (ac — eh) , Tiq = (al + ck — bh — de) / (ac — eh) . 
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While we obtain the conjugate pair propagation equations (II) from the pair propaga- 
tion equations (I) by replacing Pi,Pi + ±, g„ g i — > Viifti+ii9ii9i > that is, by replacing ratios 
rj, r i+ i, Si, St — > fj, fj + i, s i and c <-> d, h <-> I. Explicit forms are given by 

ff + rf +1 - ri(f^ + 1) - rw i+1 + r 3 f,(i + f^) + r 4 f i+1 (i + f?) = o (is) 
si + s ~f - T^sfsf + 1) - r 6 s iS ~/ - r 8 §;(i + 5 2 ) - r 7 Si(i + ~4) = o. (19) 



3. Connection between the Yang-Baxter and the Pair Propagation Equations 

Next we consider the connection between the Yang-Baxter and the pair propagation 
equations. We consider products of three i?-matrices in Eq. (Q) as matrices with row indexed 
by (3, p, p and with column indexed by a, v, a. We denote quantities in left-hand side as 
A(/3, p, p\ot, v, a) and those in right-hand side as B{(3, p, p\a, v, a). To show the relations 
A{(3, p, p\ot, v, a) = B(/3, p, p\a, v, a) is equivalent to show 

A ((3, p, p\a, v, a)v x (a)v 2 (v)v 3 (a) = Y B {(3, P, p|«> v, a)vx{a)v 2 {y)v?,{<j). (20) 
for arbitrary three vectors Vi(a)V2(v)vs(a) . Explicit forms of these equations are given by 

ri,C,<f>,a,v,cr 

= Y i? "(/ i ) p; r ?^) i? '( < / , ) C;^/?)^(^,a;^C)^i(")^2(^)^3(^) ) (21) 

rj£ ,<f> t a,is,a 

We can transform these by using the pair propagation equations as 

(left — hand side) = ^ R(p,(;r], (3)R (p,a; 4>,()vi(a)u 2 (r])u 3 ((f)) 

ri,C,<f>,ct 

= Y,R{^C\v^)A{QuMh{p) 

= Zl '(f3)z 2 \p)t 3 '(p), (22) 
(right - hand side) = ^ R " (p, P\ V, 4>)R\^, °, P) Ui' {C,)u 2 " {rj)v 3 {a) 

= Y2 R "(p,p;v,<i>)tx"W)u 2 "(vW'(<i>) 

= t/'(/?)z/( M )z 3 "(p). (23) 
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Ratios of vectors change in the following way; 



(left — hand side) : 



v 2 {-)/v 2 {+) 
{ v 3 (-)/v 3 (+) 



t'i(-)/t'i(+) ~- z' 1 (-)/z' 1 (+), 
4HA4(+) — > z' 2 (-)/z' 2 (+), 



(24) 



(III) 



-)/%(+) 



(III) 



(right — hand side) : 



: < 



)M+) 
^'(-)/4'(+) 
4'(-)M'(+) 



0? 




-)K(+) 


0i) 


of 




-)A4'(+) 


of 






*s(- 


-)/4'(+) 


(HI) 



*a(-)/*3(+), 

*"(-)/<(+), 
4'(-)/4'(+)> 
4'(-)/4'(+), 



(25) 



where we use equations (I), (II) and (III), which connect in variable X with out variable Y 
in the following forms; 

' (I) : X 2 + Y 2 - T^Y 2 + 1) - Y 2 XY + r 3 Y(l + X 2 ) + T 4 X(1 + F 2 ) = 0, 

(II) : X 2 + Y 2 - T 5 (X 2 F 2 + 1) - T 6 XY + T 7 Y(1 + X 2 ) + T 8 X(1 + Y 2 ) = 0, 

(III) : X 2 + Y 2 - T 5 (X 2 Y 2 + 1) - T 6 XF - T S Y(1 + X 2 ) - T 7 X(1 + Y 2 ) = 0. 

If the forms of Eqs. (I) and (II) are the same, Eq. fl2lD is satisfied. Conditions that the forms 
of Eqs. (I) and (II) are the same lead r, = r i+ 4, (i — 1 ~ 4). We call these conditions as the 
candidate conditions to satisfy the Yang-Baxter equations, because these conditions do not 
guarantee to satisfy the Yang-Baxter equations but mean to satisfy three vectors multiplied 
forms of the Yang-Baxter equations. 

In the 16-vertex model case, the above candidate conditions lead further restrictions on 
the Boltzmann weights, that is, we obtain two possibilities i) a + d = b + c, e = h, k = I 
and ii) e = I, k = h . Taking into account of these considerations to find the candidates, we 
consider the following exactly solvable cases, which are the more restricted cases than the 
above possibilities, i)a = c, b = d, e = h, k = I and ii)a = d, b = c, e = I, k = h in the 
next section. 



4. Exactly Solvable Cases in 16- Vertex Models 

The whole regime of 16 vertex model is not yet solved exactly, and further, not yet 
shown to be integrable. In this situation, we suspect that the 16 vertex model will be not 
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integrable in the whole regime. Then in the following, we will restrict the original model to 
more specialized cases, expecting to find integrable cases. This means that we will not cover 
the whole regime of the original model, but will cover the entire temperature range as we 
will find the whole spectrum of the transfer matrices. Furthermore, as our main poit of our 
paper is the relation between Yang-Baxter and pair propagation equations, we will examine 
some special cases of the original model in order to demonstrate the method of section 3, 
but we will not intend to analyze the whole regime completely. 

i) a = c,b = d,e = h,k = I case 

In this case, we find that the Yang-Baxter equations are satisfied in a sense that there al- 
ways exists a non-trivial set {a", , b , , e , k } for given sets {a, ,b, , e, k}, {a, ,6, , e, k }. 
By explicit calculations, apparent independent Yang-Baxter equations reduce from (32—4) = 
28 to 3 in the forms; 

e"C a - k"C p = 

(a" - b")Cfs - e'C s = (26) 
a"C 6 - b"C v = 

where C a = ak' + a k + be + b'e, Cp = ae' + a e + bk + b'k, C$ = (a — b)(a — b') + (e — 
k)(e — k ), C v = aa 1 + bb + ee + kk , = ab' + a b + ek + e k. Then for given sets 
{a, ,b, , e, k}, {a , ,b , , e , k }, there always exists a non-trivial set {a", ,b , , e , k }, 
that is, the Yang-Baxter equations are always satisfied. 
In this case, Eqs. (0), flU), ©, flO) S ive 



f (rf-l)(rf +1 -l) = 0, 

,2 



, s- 
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IK*; -i) = o, 



rf + si - T 9 ( ri 2 Si 2 + 1) - T 1Q r iSi + T u {nil + sf) - Si (l + r?)) = 0, 
^ sf + rf +1 - r 13 ( S ; 2 r i+1 2 + 1) - r 14 r iSi + T 15 (^(1 + rf +1 ) - r i+1 (l + sf)) = 0, 

where 



(27) 



T 9 = (b 2 - e 2 )/(a 2 - k 2 ), T 10 = 2(a 2 + e 2 - 6 2 - k 2 )/(a 2 - A; 2 ), T n = 2(bk - ae)/(a 2 - A; 2 ) 



l3 = (b 2 - k 2 )/(a 2 - e 2 ), T 14 = 2(a 2 + k 2 - b 2 - e 2 )/(a 2 - e 2 ), T 15 = 2{be - ak)/(a 2 - e 2 ' 



Solutions are the combinations of = Sj = ±1 and r i+ \ = s\ = ±1, (signs are independent 
for both). Then we obtain following cases; 
ia) Ti = Si = r i+ i — s i — ±1 case 

Choosing pi = p i+1 , eigenf unctions at this site become 



9i= \ 1„T!\ . g' l = (a + b±(e + k))g l , £ = 0. (28) 
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ib) ri = Si = — r i+ i = —s' i = ±1 case 

Choosing pi = p i+ i , eigenf unctions become 

* = ( ± g % ) . « = (- » ± (« - *)) ( ) ■ a = o- (») 

When we mix these cases, it is necessary to take n = rjv+i for periodicity, but otherwise we 
can mix ia) and ib) in such a way as (the number of times of mixing = — r i+1 = 1 cases) = 
(the number of times of mixing = — r i+1 = — 1 cases). General eigenvalues of transfer 
matrices give 

A = ( a + b + e + k) mi (a + b-e- k) m2 ((a - bf - (e - A;) 2 ) m3 (±l) m3 , (30) 

by using non- negative integers mi, mi and m.3, where mi + 7712 + 2m.3 = N, and mi = or 
m,2 = must be satisfied in the case 1713 = 0. 

From explicit expression of transfer matrices at small N, we obtain A = a + b±(e + k) 
for iV = 1, and A = (a + b + e + k) 2 , (a + b-e- k) 2 , ±(a - b + e - k)(a - b - e + k), for 
N=2 , which agree with the above formula. 

ii) a — d,b — c,e — l,k — h case 

In this case, the Yang-Baxter equations are not satisfied but transfer matrix commute. 
We first explain why the Yang-Baxter equations are not satisfied in this case through the 
8-vertex case, because the expression becomes rather complicated in the 16-vertex case, but 
the mechanism is the same. Then we consider this special 8-vertex case, that is, a — d, b — 
c,e = h = k = l = and the explicit Yang-Baxter equations give 

ab a + aa a = ab a + bab (31) 
abb +aab = ba b + bb b" (32) 



From Eqs.(pHD,(3"2"), we obtain (a + b')(a" — b') = , but as a + b' 7^ because Boltzmann 
weights must be positive, and as it is in general a 7^ b , the Yang-Baxter equations are not 
satisfied in general. Same mechanism happens in this 16-vertex case, and as it is in general 
a + d 7^ b + c, e + I 7^ k + h, the Yang-Baxter equations are not satisfied in this 16-vertex 
case. 

In this 16-vertex case, Eqs.(Trj), (0), (0), (0) give 

= 0, 
= 0, 
0, 

= 



-l)(rf + i-l) 

, -1)W"1) = 
l 2 -l)(^-l) = 

2 UK 1. 



(33) 
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Then we substitute solutions of Eq.(|33|) into the original pair propagation and conjugate 
pair propagation equations, and we obtain following cases ; 

iia) Ti = r i+ i = —s^ = —s i = ±1 case 

Choosing pi = p i+ i , eigenfunctions at this site becomes 

*=(=^(+))' ^ = °' 9t =(a + b T (e + k)) gi (34) 

iib) Ti = —r i+ i = Si = —s\ = ±1 case 

In this case, choosing pi = p i+1 , eigenfunctions at this site become 

*=U!m)> 9:=o ' * =( -° +6± ( e -* )) (™w) (35) 

If we mix these cases, we obtain exactly the same formula Eq . (j30l) . From the explicit 
expression of the transfer matrix of small N, we obtain A = a + b ± (e + k) for N = 1 
and A = (a + b + e + k) 2 , (a + b-e-k) 2 , ±(a — b + e — k)(a — b — e + k) for N = 2, which 
agree with this formula. 



5. Summary and Discussion 

We have clarified the connection between the Yang-Baxter and the pair propagation 
equations in the 16-vertex models. In the 16-vertex models, we find exactly solvable example 
of i) a = c,b = d,e = h,k = I case, where the Yang-Baxter equations are satisfied. 

We find another exactly solvable example of ii) a = d, b = c, e = I , k = h case. By explicit 
calculation, we can find that conditions a + d — b + c, e + l — k + h are necessary to satisfy 
the Yang-Baxter equations in this case, but these are not satisfied in general, that is, the 
Yang-Baxter equations are not necessary condition for the solvability. Though the Yang- 
Baxter equations are not satisfied, we can show that transfer matrices commute (integrable) 
for any lattice size N, which will be discussed in a separate paper. 

In the 16-vertex models, from these exactly solvable examples, integrable cases which 
satisfies the Yang-Baxter equations are rather limited cases in the whole exactly solvable 
cases. In this sense, the pair propagation equations are more fundamental, and even if the 
Yang-Baxter equations are not satisfied, if the pair propagation equations are solvable, it is 
sufficient for our purpose to find eigenvalues of tansfer matrices. 
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